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ABSTRACT 

The new extremely large telescope projects need accurate evaluation of the candidate 
sites. In this paper we present the astroclimatological comparison between the Paranal 

■ Observatory, located on the coast of the Atacama Desert (Chile) , and the Observatorio 

del Roque de Los Muchachos (ORM), located in La Palma (Canary Islands). We 

Q . apply a statistical analysis using long term databases from Paranal and Carlsberg 

' Meridian Telescope (CAMC) weather stations. The monthly, seasonal and annual 

C/2 | averages of the main synoptical parameters in the two sites are computed. We compare 

the long term trends in order to understand the main differences between the two sites. 
Significant differences between the two analyzed sites have been found. Temperature 
have increasing trends in both observatories with somewhat higher evidence at ORM. 

■ Seasonal variations of pressure at Paranal have been highly decreasing since 1989 and 
| we do not see the same phenomenon at ORM. The two sites are dominated by high 

pressure. In cold seasons RH is lower than 60% at CAMC and 15% at Paranal. In 
warm seasons RH is lower than 40% at CAMC and 20% at Paranal. The analysis of 
| the dew point has shown better conditions at Paranal with respect to CAMC in winter, 

\Q . autumn and spring before 2001, while the two sites are becoming similar afterwards. 

' Winds at ORM are subject to pronunced local variations. 

| Key words: Site testing - Atmospheric effects - Methods: data analysis. 

2 i 1 INTRODUCTION cillating weather variation due the presen ce of El Nino and 

, , „ . , La Nina events (Sarazin(author?) 1 19971 ). the ORM is in- 

This study follows a series of papers aimed to compare n , , . , , ... 

ffuenced by a semipermanent Azores high pressure system 
two oi the most important sites candidates to host the fu- , ., . . „ , , , .... . ,. . ,, 

and it is influenced by the almost periodicaf variation of the 
turc extremely large ground based telescopes. The analysis AT ... „ .„ .. ,„,„,,,„ , . , 

. , . . ,.. North Atlantic Oscillation (NAO) (Wanner et al. (author?) 

or the astronomical and meteorological parameters m dif- „ , , .. „,. L„„J T , ,. . . ; ox 

° . 2001; Graham(author :) 120051 ; Lombardi et al. (author;) 



ferent sites are conducted, since today, using different ap 
proach. In this work we want to apply same methods to 
compare two different and important sites for the develop- 
ment of the astronomical observations: the site of Paranal 
Observatory located in the Atacama Desert (Chile, south- 



12009) . 

More than 20 years of meteorological data have been col- 
lected at ORM using the Carlsberg Automatic Meridian 
Circle (CAMC) meteorological station and a detailed anal- 
ysis can be found in Lombardi et al. (200q) (hereafter Pa- 
MucilZf^RMUoc PCT T > and L ° mbardi 6t aL » (hereafter Paper II). Pa- 

IVlUCIlclCIlOS I vylXlvl I lOCdieU <Al LitX r^dllllcL I K^dulldL V ISlcLIlQS, x . . 1-^,1 , ■ 1 , 

v ' v J per 1 presents a complete analysis of the verticaf tempera- 

northern hemisphere). The two opposite locations have dif- ,. .mi . n T • 1 /-< 1-1 ^^ T ^\ i 

, , ture gradients at lelcscopio INazionale Cahleo (1ING) and 

ferent observing characteristics not only because it is pos- ,, . , ,. . , . T , i t-< 

their correlation with the astronomical seeing. Instead, Pa- 
sible to see two different sky regions, but also because the TT , . r , .. , , . , 

J . per 11 reports an analysis ot the correlation between wind 

oceanic currents impress different climatic regimes to the , . , ,, , 

and astronomical parameters as well as the overall long term 
two sites. While Paranal is located m a desert area, 12 km , . „t-.hi 

weather conditions at ORM. Differences m the ORM micro- 
far from the Pacific Ocean, and it is subjected to the os- ,. , , , . , . , . , 

J climate have been demonstrated in a detaifed comparison 

between synoptical parameters taken at three different lo- 
cations in the observatory on a 1000 m spatial scale. ORM 

* E-mail: glombard@eso.org 
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Table 1. Geographical positions of Paranal and CAMC. 





Latit 


ode 




Long 


itudc 


Height [m a.s.l.] 


Paranal 


24° 37' 


31" 


S 


70° 24' 


10" W 


2636 W 


CAMC 


28° 45' 


36" 


N 


17° 52' 


57" W 


2326l b ' 



l a l Platform altitude. 
' b l Dome floor. 



Table 2. Available databases for the observatories. The data cov- 
erage with respect to the Total is also reported. 





TNG 


CAMC 


Paranal 


Data rate 


30 sec 


5 min 


20 min 


Begin 


March 1998 


May 1984 


January 1985 


End 


December 2007 


March 2005 


December 2006 


Total 


~10 yr 


~21 yr 


~22 yr 


Data coverage 


T 


87% 


87% 


80% 


P 


87% 


86% 


71% 


RH 


85% 


85% 


80% 


Tr>p 


85% 


85% 


80% 


w dir 


87% 


85% 


75% 


VJ B p 


87% 


86% 


80% 



shows to be dominated by high pressure, and characterized 
by an averaged r elative humidity lower than 50%. Finally, in 
Lombardi et al. (|2008al ) we have analysed in detail the prop- 
erties of the dust concentration on La Palma from ground 
based measurements and estimated the aerosol extinction in 
B, V and / on the site. 

Extensive site testing campaigns have been conducted on the 
top of the Paranal Observatory since years. Thanks to the 
excellent results, the site was chosen to host the four Very 
Large Telescopes (VLT) by the European Southern Obser- 
vatory (ESO). 

This site, like La Silla, the other Chilean ESO site, have been 
very deeply analyzed by ESO teams. Now ESO telescopes 
are considered the touchstones and their characteristics are 
used to be compared with the other sites and the other tele- 
scopes. 

The present paper is organized as follows: 

• Section 2: describes the Paranal Astronomical Site 
Monitor and the data reduction; 

• Section 3: compares the temperature, pressure, rela- 
tive humidity and correlate each parameter with the South- 
ern Oscillation Index (at Paranal) and NAO (at ORM); 

• Section 4: shows the dew point comparison at the two 
sites; 

• Section 5: analysis of wind direction and wind speed; 

• Section 6: summarization of the final conclusions. 



2 DATA REDUCTION 

Table [1] that lists positions and heights of the telescopes at 
observatories. Paranal is almost 300 m higher than CAMC 
and about 4° closer to the Equator. 

Table [2] reports the databases characteristics. The Paranal 
Astronomical Site Monitor is located in the north area of 



the Paranal Observatory platform and hosts several instru- 
ments used in the characterization of the site. In particular 
in this paper we make use of the meteorological data from 
the Vaisala tower. 

The Vaisala tower is a robust steel structure having a to- 
tal height of 30 m. All the data are regularly collected since 
January 1st, 1985 and have to be intended as 20 minutes av- 
erages. The external air temperature (T) is simultaneously 
measured at 2 and 30 m above the ground (same height of 
the VLT's domes) with an accuracy of ±0.2° C (Sandrock et 
al.(author?) ll999l l. Air pressure (P) and relative humidity 
(RH) are measured at 2 m above the ground with an accu- 
racy of ±0.1 hPa and ±1% respectively. 
Wind speed (w sp , in [m s -1 ]) and wind direction measure- 
ments are taken at 10 m above the ground with a precision 
of ±0.5 m s _1 for the wind speed and ±3° for the wind di- 
rection. Following the same statistical procedures of Paper 
I and Paper II we have computed the hourly averages and 
then the monthly averages starting from the T, P, RH and 
w S p raw data series. 

A particular care is used to minimize any effect due to biases 
in case of lacking of data that typically occurred in winter 
time. For each missing month value we take into account 
the average obtained from the two corresponding months in 
other years in which the values of the months before and 
after the absent one are similar. For example, if the lack- 
ing month is September 2002, we look for the two Augusts 
and Octobers in the other years having similar mean values 
of August and October 2002. The accepted September 2002 
value is the average of the Septembers corresponding to the 
chosen Augusts and Octobers. This is the main reason why 
we decided to use monthly averages as an intermediate step 
in the calculation of the annual averages. Finally, we com- 
pute the annual averages of T from the monthly ones for the 
three telescopes (see Paper I). 

We provide to report in tables the annual values and the 
rms of the annual averages of the most important meteo pa- 
rameters. 

Wind direction statistics is evaluated by calculating the an- 
nual percentage of hours in which the wind comes from each 
direction D. The wind rose has been divided into 8 mean 
directions (N, NE, E, SE, S, SW, W, NW) and the per- 
centages of hours are calculated into intervals defined as 
[_D- 22.5°, Z)± 22.5° [. 

CAMC data from Paper I and Paper II are used in order 
to compare ORM with Paranal. For what concerns the tem- 
perature trend at ORM we decided to take into account also 
data recorded at TNG because they extent to year 2007 (3 
years more than CAMC). Telescopes at ORM are located 
on a space baseline of about 1000 m (see Paper I). 
The TNG meteo tower is a robust steel structure with a total 
height of 15 m. The tower is located about 100 m far from 
TNG building. The data are regularly sent from the tower 
to TNG annex building by means of an optic fiber link since 
27 March 1998. The data sampling rate is 10 seconds, while 
data storage is done every 30 seconds. 

The CAMC carried out regular meteorological observations 
in the period 13 May 1984 to 31 March 2005 and the records 
are more or less continuous in that period. For the years 
1984, 1985 and 1986 meteorological readings are only avail- 
able at 30 minute intervals. From January 1987 readings 
were made at 5 minute intervals throughout the day and 
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Figure 1. Annual temperatures at Paranal (circles), CAMC 
(squares) and TNG (filled circles). The solid and the dotted lines 
indicate the linear fit of the data for Paranal and CAMC respec- 
tively. 

night regardless of whether observing was in progress. Be- 
ginning in December 1994, all readings were made at 20 
seconds intervals and then averaged over 5 minutefl 
Both Paranal and ORM are located well above the inversion 
layer, in fact the altitude of the invers ion lay er at Paranal is 
about 1000 m as reported by Sarazin (1994), while at ORM 
it occurs in the rang e betw een 800 m and 1200 m (Mclnnes 
& Walker(author?) Il974l ). 

3 ANALYSIS OF T, P, SOI AND RH 
3.1 Temperature 

Table|3]reports the computed mean annual temperatures for 
Paranal, CAMC and TNG. Values in parenthesis correspond 
to the rms of the annual averages. 

The Paranal temperature is taken at two different levels, 
at 2 and 30 m above the ground (Tb and T30 respectively). 
The vertical variation of temperature with the altitude (wet 
and dry adiabatic lapse rate) i s betw een 6.0 and 10.0° C/km 
(Kittel & Kroemer(author?) Il980l ). so we expect a ~0.2- 
0.3°C difference in temperature between the two sensors at 
2 and 30 m. As shown in Table [3] the mean difference of 
the temperature taken at the two heights is 0.2° C, that is 
comparable with the accuracy of the sensor. For this reason 
we decide to use the 22 years long database at 2 m as repre- 
sentative of the temperature at Paranal. For completeness, 
Table [3] reports also the annual T of CAMC measured at 
10.5 m above ground and TNG measured at 10.0 m above 
ground (see Paper I). 

Figure [1] shows the plot of the Tb annual mean values re- 



1 http://www.ast.cam.ac.uk 



Figure 2. Thermal excursion distribution at Paranal (solid line) 
and CAMC (dotted line) after 1998. The two sites are almost 
equivalent. 



ported in Table [3] as well as the annual temperature mea- 
sured at CAMC and TNG. The data show an offset of 
about 4°C between Paranal and ORM. Furthemore the 
rms of the annual averages are about 3 times higher at 
ORM with respect to Paranal. This suggests an higher tem- 
perature variations during the years at ORM. Trends at 
Paranal and CAMC show a positive slope during the years, 
while TNG data have a flat trend. The best linear fit of 
CAMC data gives an increase of the temperatures of about 
(1.0 ± 0.3)°C/10yrs, while the slope computed for Paranal 
gives a value of (0.4 ± 0.1)°C/10yrs. This is not surprising 
considering that the global warming and the glacier retreat 
are less pronounced in South America. To facilitate the com- 
parison, it is also drawn the total average for Paranal and 
CAMC. 

But it is interesting to note that after 1998 data at TNG 
seem to indicate a flattening of the increasing temperature 
trend. This rises the question if we are in presence of a global 
warming or it is a typical temperature oscillation through 
decades and/or a regional effect. 

The 22 years long baseline of Paranal is characterized by 
an increasing trend, the linear regression of Paranal data in 
the period 1993-2000 shows a slope of (-0.12 ± 0.05) while 
data in the period 2000-2006 have an opposite slope slightly 
steeper (0.19 ± 0.07). This suggests a possible correlation 
with the occurrence of wide-scale climatological events such 
as El Nino and La Nina phenomenons. In fact, a preliminary 
check has shown the presence of strong La Nina episodes in 
1999 and 2000 that are linked to the coldest annual temper- 
atures at Paranal (see Section T3.4p . 

This confirms the strong link between air and ocean temper- 
ature that may influences the high level of the atmosphere. 
The year 1989 appears to be the coldest for both sites in 
the last 20 years. The year 2001 is the warmest for CAMC 
while the 2003 and 2006 are the warmest for Paranal. 
The typical 3—4 years oscillation of the temperatures at 
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Table 3. Mean annual temperatures at Paranal, CAMC and TNG [°C]. Values in parenthesis correspond to the rms of the annual 
averages. 



Year 


1985 


1986 


1987 


1988 


1989 


1990 


1991 


1992 


Paranal T2 
Paranal T30 
CAMC Tio 
TNG T 10 


12.2(1.4) 
8.8(4.4) 


12.7(1.1) 
8.9(3.9) 


9.1(3.1) 


12.8(2.0) 
7.4(5.4) 


11.9(1.6) 

5.2(4.3) 


12.2(1.7) 
8.8(5.1) 


8.7(4.8) 


7.9(4.2) 


Year 


1993 


1994 


1995 


1996 


1997 


1998 


1999 


2000 


Paranal T2 
Paranal T30 
CAMC Tio 
TNG Tio 


12.8(1.4) 
7.0(5.2) 


13.2(1.3) 
9.8(4.9) 


13.1(1.3) 
9.5(4.3) 


12.5(1.2) 
8.6(5.1) 


13.1(1.6) 

8.9(4.0) 
10.1(5.1) 


12.9(2.1) 

10.0(5.0) 
9.6(5.3) 


12.3(1.4) 

9.3(5.2) 
9.9(5.3) 


12.0(1.7) 
11.9(1.7) 
9.6(5.1) 


Year 


2001 


2002 


2003 


2004 


2005 


2006 


2007 


Average 


Paranal T2 
Paranal T30 
CAMC Tio 
TNG Tio 


12.8(1.3) 
12.6(1.3) 
10.1(4.3) 
10.7(4.1) 


12.8(1.6) 
12.6(1.7) 
9.6(4.3) 
9.7(4.5) 


13.6(1.3) 
13.5(1.3) 
9.8(4.8) 
9.7(4.8) 


12.9(1.5) 
12.7(1.4) 
9.0(4.9) 
8.9(5.0) 


13.0(1.3) 
12.6(1.4) 

9.5(4.1) 


13.6(0.8) 
13.3(0.9) 

10.0(4.7) 


9.6(5.1) 


12.8 ± 0.5 
8.8 ± 1.2 



Table 4. Mean annual air pressure at Paranal and CAMC [hPa]. Values in parenthesis correspond to the rms of the annual averages. 



Year 


1985 


1986 


1987 


1988 


1989 


1990 


Paranal 
CAMC 


774.3(2.6) 


775.5(2.3) 


774.1(1.8) 


773.8(2.2) 


743.3(1.2) 
773.6(2.3) 


743.3(0.9) 
774.2(2.4) 


Year 


1991 


1992 


1993 


1994 


1995 


1996 


Paranal 
CAMC 


774.4(2.5) 


774.5(2.2) 


743.1(0.7) 
773.9(2.7) 


743.5(0.5) 
775.1(1.8) 


743.7(0.7) 
774.6(2.2) 


743.3(0.6) 
775.5(3.8) 


Year 


1997 


1998 


1999 


2000 


2001 


2002 


Paranal 
CAMC 


743.3(0.7) 
775.4(2.4) 


743.5(0.5) 
775.9(2.2) 


743.5(0.4) 
775.4(2.2) 


743.2(0.6) 
775.4(2.4) 


743.4(0.4) 
776.2(2.5) 


743.5(0.3) 
776.1(2.4) 


Year 


2003 


2004 


2005 


2006 




Average 


Paranal 
CAMC 


743.0(0.4) 
776.1(2.4) 


743.6(0.4) 
776.0(2.7) 


743.8(0.4) 


743.8(0.4) 




743.5 ± 0.2 
775.0 ± 1.0 



CAMC pointed in Paper I is not so clear for Paranal. 
An interesting results comes from the analysis of the daily 
thermal excursion after 1998 in the two sites, in fact looking 
at Figure [2] we do see a very similar distribution for both 
sites with a peak at bin number 18, corresponing to a ther- 
mal excursion of about 4.8° C. 



3.2 Air pressure 

Figure [3] shows the annual mean pressure P calculated at 
Paranal and CAMC. Table [4] reports the values plotted in 
the figure. The average pressure at Paranal during the entire 
baseline is 743.5 hPa, while the average pressure at CAMC is 
775.0 hPa. The plot shows two different behaviors at Paranal 
with a changing of slope in the year 2000. Splitting the fit in 
two subsamples we see that the best fit of the points before 
the 2000 is almost flat (±0.1 hPa) with a mean value of 743.4 
hPa. The best fit after the year 2000 is characterized by a 
positive slope of about (0.1±0.1) hPa and the mean pressure 
is 743.5 hPa. This difference is comparable to the accuracy 
of the measuring instruments. A more detailed analysis is 
needed to check if this is a long term increasing trend of 
pressure. The data obtained from CAMC show a similar in- 



creasing trend during the years (~1.5 hPa, see Paper II). 
The minimum pressure measured at Paranal was 728.0 hPa 
in September 1989, while the highest one was in April 1998 
with a value of 754.0 hPa. Instead, CAMC minimum pres- 
sure was 747.7 hPa in September 2004, while the maximum 
was 785.8 hPa in July 2001. 

Two interesting plots are shown in Figures [4] and [5] where 
the monthly P averages in the two sites are plotted. There 
is a clear decreasing of the dispersion of the monthly P dur- 
ing the years at Paranal (Figure [4] top) as confirmed by the 
analysis of the standard deviation of the annual averages 
(Figure [4] bottom). In fact at Paranal it decreased of about 
70% between 1990 and 2006, while the same data computed 
at CAMC do not show the same phenomenon (Figure [SJ). 
The same analysis applied to monthly temperatures shows 
a lower seasonal dispersion at Paranal (~8°C) with respect 
to CAMC (~19°C). Temperature dispersion remained al- 
most constant up to now in both sites. 

To check if the effect of the yearly variation of the pressure 
at Paranal can be explained taking into account wide-scale 
climatological episodes, we have checked a possible correla- 
tion with the Souther Oscillation Index. This analysis will 
be discussed in Section 13.41 
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Table 5. Seasons definition at Paranal and CAMC. 





Paranal 


CAMC 


WINTER 


J uly- August- Sept ember 


January-February-March 


SPRING 


Octobcr-Novcmbcr-Dcccmbcr 


April-May- June 


SUMMER 


January- February- March 


July- August- September 


AUTUMN 


April-May- June 


Octobcr-Novcmbcr-Dcccmbcr 



1995 2000 

YEAR 




1995 
YEAR 



g 744 - 

W 

m 

g 743 - 
5 740 - 
738 — 



H 1 1 1 1 1 1 1 1 H 



H 1 1 1 1 H 



1990 1995 2000 

YEAR 



Figure 3. Top: annual P at Paranal Observatory. Bottom: an- 
nual P at CAMC. The averages are represented by the dotted 
lines, while the linear fits by the solid lines. 



Figure 4. Top: Monthly P at Paranal. Bottom: Standard devi- 
ation of the monthly P at Paranal. 



Following Paper II, we have calculated the theoretical pres- 
sure for Paranal using the barometric correction that de- 
pends on sites scale height H in the barometric law. The 
knowledge of the average pressure and the atmospheric scale 
height are important parameters to compute the local con- 
tribution of the Rayleigh scattering to th e tota l astronomical 
extinction (Hayes & Latham(author?) Il975h . 
At CAMC we have found a typical scale height Hcamc = 
8325 m and a theoretical pressure of 766.0 hPa in the pe- 
riod between 1998-2004, while the data give 775.9 hPa on 
average, confirming that ORM is dominated by high pres- 
sure (Paper II, Section 5.1). Since the barometric correc- 
tion is a function of the mean temperature (Ti ayer ) of the 
layer in which the correction is applied, we have taken into 
account as closest weather station the city of Antofagasta, 
close to the Pacific Coast, few meters above sea level. Be- 
tween 1998 and 2004 the mean temperature of Antofagasta 
was 16.7 ± O.4°C0, while the mean temperature computed 
at Paranal in the same years was 12.8 ± 0.5°C. These two 
mean temperatures (Paranal and Antofagasta) allow us to 
calculate the mean temperature of the layer between the sea 
level and the Paranal, using the weighted average of the two 
measurements. We obtain {Ti ay er) = 15.2 ± 0.5°C. 
The standard atmospheric model from NASA's Glenn Re- 



1990 1995 2000 

YEAR 



Figure 5. Top: Monthly P at CAMC. Bottom: Standard devia- 
tion of the monthly P at CAMC. 



2 See http://www.tutiempo.net 
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Figure 6. Averages of monthly P (top) and RH (bottom) at 
Paranal. 



Figure 7. Averages of monthly P (top) and RH (bottom) at 
CAMC. 



search Center (GRC) Web sit^fl permits us to calculate 
the theoretical pressure at the altitude of Paranal for a 
(Tiayer) = 15.0° C which gives a resulting theoretical pres- 
sure of 734.4 hPa, and a theoretical scale height H grc ~ 
8238 m. Instead, using standard tables in Allen (2000), 
for (Ti ayer ) = 15.0° C, we find a theoretical scale height 
H Alien ~ 8430 m, in good agreement with Hgrc- The mean 
value between Hgrc and Hah en gives US H Paranal — 8334 
m which corresponds to a theoretical pressure for Paranal of 
Pparanai = 738.3 hPa. This result is lower with respect to 
the empirical pressures reported in Table [4] so we can con- 
firm that also Paranal is dominated by high pressure. The 
Paranal scale height is a bit higher than at ORM (8334 m vs 
8325 m) as expected having a lower latitude site and higher 
average temperature than ORM. 



3.3 Relative humidity 

The relative humidity (RH) and the dew point are two im- 
portant parameters for the astronomical instrumentation, 
because they set the occurrence of moist and water con- 
densation on the coldest part of the telescope and of the 
instruments. In particular these parameters may affect the 
upper surface of the main mirror and the pipes of the cool- 
ing system. 

We separated the annual RH in 4 seasons as defined in Table 
Note that Paranal and ORM are in opposite hemispheres, 
so the definitions are inverted for the two observatories. 
To facilitate the discussion we consider as cold season au- 
tumn and winter, and warm season spring and summer. Ta- 
bles [()] and [7] report the computed annual RH in cold and 
warm seasons at Paranal and CAMC. Typically in cold sea- 
sons RH is lower than 60% at CAMC and 15% at Paranal. 
In warm seasons RH is lower than 40% at CAMC and 20% 



See http: / /www.grc. nasa.gov 



at Paranal. Figures [6] and [7] show the seasonal trend of P 
(top) and RH (bottom) in both sites considered as the aver- 
age of the same month through the considered years of the 
databases. The variation of P at Paranal is < 1 hPa, very 
lower with respect to CAMC (~6 hPa). Variations of pres- 
sure in a short time scale (few hours) can induce weather 
instabilities, while in a long time scale they increase the dif- 
ferences in the weather in different seasons. 
We clearly see an anti-correlation between P and RH in 
summer at CAMC, while at Paranal this effect is not ob- 
vious. Because of the effect of the Bolivian Winter, driving 
equatorial humid air from amazonian basin along the Andes, 
the RH at Paranal is higher in warm seasons (~20%) with 
respect to cold ones (~12%). Higher RH appears clearly in 
wintertime than in summertime at CAMC. 
A standard requirement for the use of telescopes is a RH 
value < 80% or < 85%. For both sites we have calculated 
the number of nights in which RH has been higher than 
the mentioned limits for more than 50% of the duration of 
the night. Only nights which duration has been ^ 6 hours 
have been used in the calculation. Results are reported in 
Table [8] and denote a significant difference between Paranal 
and CAMC, the first appearing almost immune to high RH 
events. The number of nights at CAMC is the same for the 
two imposed limits. 



3.4 The Southern Oscillation Index (SOI) 

In Paper I we demonstrated the correlation between the 
North Atlantic Oscillation Index and the annual mean tem- 
peratures at CAMC (c.l. ~0.9). The NAO is the dominant 
mode of atmospheric circulation in the North Atlantic re- 
gion and is generally defined as the difference in pressure 
between the Azores high pressu re and the Icelandic low 
press ure (Wanner et al. (author?) 2001; Graham(author?) 
120051 : Paper I Section 4). We pointed that the action of a 
positive NAO Index is as a brake for the increase in tern- 
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Table 6. Mean annual RH in cold seasons at Paranal and CAMC [%]. Values in parenthesis correspond to the rms of the annual averages. 



averages. 



Year 


1985 


1986 


1987 


1988 


1989 


1990 


1991 


1992 


Paranal 
CAMC 


9.1(2.4) 
40.3(7.0) 


9.7(4.6) 
43.8(4.1) 


48.2(2.5) 


9.2(5.1) 
52.2(7.0) 


12.0(5.5) 
49.2(3.6) 


9.6(2.1) 
56.5(14.7) 


57.4(10.2) 


52.7(4.7) 


Year 


1993 


1994 


1995 


1996 


1997 


1998 


1999 


2000 


Paranal 
CAMC 


10.8(3.8) 
47.8(9.3) 


11.0(4.2) 
40.4(8.8) 


10.5(3.1) 
40.6(8.6) 


10.7(3.0) 
53.1(12.7) 


14.8(3.9) 
44.7(12.4) 


12.6(4.3) 
30.9(6.1) 


10.9(2.0) 
12.5(7.7) 


14.9(4.7) 
53.0(13.6) 


Year 


2001 


2002 


2003 


2004 


2005 


2006 




Average 


Paranal 
CAMC 


12.5(4.3) 
33.7(14.7) 


14.5(6.3) 
56.4(8.9) 


11.4(3.0) 
36.3(10.8) 


12.4(3.0) 
40.8(12.8) 


11.5(1.3) 


12.4(3.6) 




11.6 ± 1.8 
44.5 ± 10.1 


n annual RH in warm 


seasons at Paranal and CAMC [%]. Values 


in parenthesis correspond to the rms 


Year 


1985 


1986 


1987 


1988 


1989 


1990 


1991 


1992 


Paranal 
CAMC 


43.1(3.5) 


16.9(8.3) 
42.0(7.3) 


15.8(7.0) 
43.6(4.6) 


13.0(4.6) 
30.8(7.8) 


16.9(12.9) 
35.9(10.1) 


11.5(4.4) 
42.5(14.3) 


16.3(6.8) 
29.2(10.3) 


26.1(12.0) 


Year 


1993 


1994 


1995 


1996 


1997 


1998 


1999 


2000 


Paranal 
CAMC 


21.0(14.1) 


15.8(5.7) 
16.7(11.1) 


16.2(6.3) 
26.5(10.8) 


15.0(7.3) 
15.0(6.2) 


20.0(12.2) 
15.6(15.0) 


16.4(7.6) 
6.5(5.1) 


21.1(12.2) 
19.7(10.2) 


25.5(14.4) 
24.3(13.6) 


Year 


2001 


2002 


2003 


2004 


2005 


2006 




Average 


Paranal 
CAMC 


22.8(15.4) 
21.9(9.1) 


19.8(11.1) 
22.0(15.5) 


15.3(8.2) 
24.5(3.6) 


17.0(8.6) 
32.2(9.8) 


18.0(6.0) 






17.4 ± 3.4 
27.0 ± 10.3 



Table 8. Paranal and CAMC: annual number of nights in which 
RH has been higher than 80% and 85% for more than 50% of the 
duration of the night. Only nights which duration has been ^ 6 
hours have been used in the calculation. For each year, the total 
number of nights having duration 6 is also reported. 



Table 9. Strongest El Nino and La Nina episodes between 1985 
and 2006. 



RH > 80% 

Year Paranal CAMC 



RH > 85% Nights 
Paranal CAMC Paranal CAMC 



1984 




2 




2 




155 


1985 


1 


8 





8 


324 


253 


1986 





8 





8 


311 


285 


1987 





22 





22 


108 


332 


1988 





51 





51 


258 


34!) 


1989 


1 


68 





68 


310 


362 


1990 





73 





73 


329 


354 


1991 





53 





53 


177 


358 


1992 




58 




58 




366 


1993 





58 





58 


300 


363 


1994 





41 





41 


328 


355 


1995 





51 





51 


333 


355 


1996 





81 





81 


348 


354 


1997 


2 


52 





52 


336 


304 


1998 





36 





36 


348 


359 


1999 


2 


48 


2 


4K 


360 


354 


2000 


4 


49 


2 


49 


365 


333 


2001 










(id 


364 


353 


2002 


4 


31 


3 


31 


362 


274 


2003 





27 





27 


352 


346 


2004 





35 





35 


365 


334 


2005 


3 




3 




233 





peratures, and an accelerator for the decrease. Conversely, a 
negative NAO Index acts in the opposite manner (see Paper 
I, Section 4). 

In the case of Paranal we have analysed the Southern Oscil- 



El 


Nino 


La Nina 




1987 


1988 


1992- 


■1994 


1996 




1997 


1999-2000 



lation Index (SOI). The SOI is defined as the difference in 
air pressure occurring between the western and eastern trop- 
ical Pacific (Tahiti and Darwin, Australia). Changes in the 
SOI correspond also with changes in temperatures across the 
eastern tropical Pacific. The negative phase of the SOI repre- 
sents below-normal air pressure at Tahiti and above-normal 
air p ressure at Darwin (Halpert & Ropelewski(author?) 
Il992h . Prolonged periods of negative SOI values reflect ab- 
normally warm ocean waters across the eastern tropical 
Pacific, typical of El Nino episodes, while periods of pos- 
itive SOI values coincide with abnormally cold ocean wa- 
ters across the eastern tropical Pa cific, typical of La Nina 
episodefl (Higgins et al. (author?) 120011 ). 
Because of the strong influence of the SOI on the meteoro- 
logical conditions in the southern hemisphere, and on ob- 
serving conditions at Paranal and La Silla, it is important 
to investigate possible correlations between the SOI and the 
temperature and pressure at Paranal. For this reason we 
calculated the annual averages of the SOI from the monthly 
averages retrieved from the National Prediction Center Web 



See http://www.cpc.ncep.noaa.gov 



8 G. Lombardi et al. 



1 I I I I 1 , I I I 1 

^ * A 

A 


1 , i i i | 

A A 
k A A A A 
A A 

' 1 1 1 ' 1 1 1 \ 


■ — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 








A 

A I 
A A A A 

A 


i i i i i i i i ; 

A 

A A 

> A A A A 
A 



1985 1990 1995 2000 2005 

YEAR 



Figure 8. Southern Oscillation Index (middle) vs. annual tem- 
peratures (top) and annual air pressure (bottom) at Paranal. The 
error-bars in the middle plot correspond to the rms of the SOI 
annual means. 



Time. Only clear nights have been taken into account. The 
declared percentage of clear nights at ORM is 84% (Munoz- 
Tunon et al. (author ?) 120071 ). while at Paranal it is 85% 
(http:/ /www.eso.org). 

The analysis concerns the annual percentage of time in 
which AT = T — Top < X, where X corresponds to a 
variable upper limit for AT (X = 1 and 5°C). This statis- 
tics is crucial for the knowledge of the amount of time in 
which a danger of condensation on the telescopes hardware 
may occur in two extreme AT. 

In Figure 9 we report the computed frequencies in the two 
sites. In winter, spring and autumn we see a significant dif- 
ference between Paranal and CAMC until 2000. After such 
year the two sites seem to became more similar. 
Paranal shows negligible percentages for all X limits in 
spring, while it never goes above ~4% for X — 5°C in winter 
(2002) and autumn (2002 and 2004). CAMC has frequency 
often > 10% when X = 5°C. After 2001 CAMC shows sim- 
ilar percentages with respect to Paranal when X = 1°C. 
This is due to the very little increasing of the percentages 
at Paranal after 2001. 

In summer the two sites are very similar maintaining per- 
centages never higher than ~3% at Paranal and ~4% at 
CAMC when X = 5°C, with the exception of the year 1990 
for CAMC. 



site. 

Figure [8] (middle) shows the SOI as a function of the year 
where the error bars correspond to the rms of the SOI annual 
means. Each minimum in the figure corresponds to warm 
episode (El Nino), while each maximum to cold episode (La 
Nina). As shown in the figure and also reported in Table |9j 
strong El Nino occurred in the years 1987, 1992-1994 and 
1997, while strong La Nifia occurred in 1988, 1996 and 1999- 
2000. 

The top of Figure [8] shows the T annual averages at Paranal. 
We see an increase in the annual temperatures in connection 
with El Nino, vice versa minimum temperatures occur with 
La Nina episodes (correlation is ~0.7). Local minimums in 
the P annual trends at Paranal (Figure[S] bottom) are linked 
to the presence of La Nina events, while maximums with El 
Nino, but in this case poor correlations is found between the 
yearly trends at Paranal. 



4 ANALYSIS OF THE DEW POINT 

The dew point temperature (Top) is the critical tempera- 
ture at which condensation occurs. When dew point tem- 
perature and air temperature T are equal the air said to 
be saturated and condensation appears if the air cools. It 
is clear that the knowledge of the dew point temperature is 
crucial for the maintenance of the optics of the telescope to 
avoid condensation, in particular if the the instrumentation 
is maintained at temperatures few degrees lower than the air 
temperature. Condensation can be reached mostly in clear 
nights, when the earth cools rapidly, therefore it is important 
to know the percentage of time in which condensation may 
occur to have the best performances from ground based tele- 
scopes. To this aim we have computed the dew point tem- 
perature using nighttime data for both Paranal and CAMC. 
We have defined as nighttime the range 20:00—6:00 hr Local 



5 ANALYSIS OF THE WIND 
5.1 Wind direction 

In the analysis of the wind direction at Paranal we made 
use of a 9 years database between 1998 (VLT first light) 
and 2006. In Paper II the time range used in the computa- 
tion of the wind statistics at CAMC was between 1998 and 
2004. We consider only nighttime data, between 22:00 and 
4:00 hr Local Time. 

In Paper II we have demonstrated that there are signifi- 
cant differences in the wind direction and wind speed ranges 
among the ORM. In such paper we considered the locations 
of the CAMC, the Telescopio Nazionale Galileo (TNG) and 
the Nordic Optical Telescope (NOT) . Because of the differ- 
ences noticed in the three sites, none of them can be consid- 
ered as fully representative of the ORM in wind analysis, so 
in this case we need to consider them together. In Figure [TO] 
we show the wind roses computed for Paranal and CAMC, 
TNG and NOT (see also Paper II, Figure 2). 
At CAMC there is no evidence of a prevailing direction. 
In Paper II we pointed out that northern winds seem to 
oscillate with a period of 10 years, while winds from the 
north-west show a similar oscillation in the opposite phase. 
TNG shows a dominant north-east mode, while NOT has 
two dominant wind d irecti ons (west and east). 
Munoz-Tunon et al. (| 19981 ) demonstrated the homogeneity 
of the image quality among the ORM. Considering the dif- 
ference between nighttime wind roses at different locations 
in the observatory that we have mentioned above, and con- 
sidering also the complex orography of the ORM, should be 
interesting to investigate if any difference exists in the lo- 
cal Surface Layer seeing among the observatory. This is an 
important issue because local conditions have been demon- 
strated to be crucial for image quality at Paranal (Lombardi 
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Figure 9. Nighttime annual percentage of time in which AT = T — Top < 1 and 5°C at Paranal and CAMC in winter, spring, summer 
and autumn. 



et al.(author?) l2008bl and Sarazin et al.(author?) 120081 ). 

At Paranal a dominant wind blowing from north and north- 
east appears during the night. Table [TOl reports the com- 
puted wind direction frequencies at Paranal between 1998 
and 2006. The values in the table are plotted in Figures [TT] 
and[l2]where we include the wind frequencies since 1985 dis- 
tinguishing two different epochs before and after the VLT 
first light occurred in 1998. Wind from north shows a clear 
decreasing trend through the years togheter with an increas- 
ing of the wind direction from south-east. This result is in 
agreement with Sarazin (2004) that in a previous analysis 
show a progressively replacement of the north-westerly wind 
by a north-easterly wind. To better investigate this behav- 
ior we have computed the yearly evolution of the nighttime 
frequencies of the wind in each direction. There is a strong 



Table 10. Nighttime wind direction frequencies at Paranal from 
1998 to 2006. 



Year 


N 


NE 


E 


SE 


S 


sw 


W 


NW 


1998 


48.6 


23.6 


6.5 


4.3 


2.0 


0.6 


1.1 


13.3 


1999 


46.7 


23.3 


6.2 


8.0 


3.7 


0.4 


1.2 


10.5 


2000 


46.8 


21.1 


5.8 


8.9 


4.7 


1.0 


1.7 


10.0 


2001 


45.8 


19.6 


7.8 


9.8 


3.4 


1.1 


2.0 


10.5 


2002 


46.4 


21.3 


6.3 


7.7 


3.3 


1.4 


2.7 


10.9 


2003 


38.5 


19.8 


7.2 


14.6 


4.(1 


1.9 


2.0 


12.0 


2004 


42.8 


23.3 


8.7 


10.1 


3.0 


1.8 


1.5 


8.8 


2005 


38.6 


19.3 


5.6 


12.3 


8.0 


3.4 


1.1 


11.7 


2006 


35.0 


26.7 


7.8 


12.3 


6.5 


2.0 


1.5 


8.2 


Tot. 


43.8 


22.0 


6.8 


9.6 


4.2 


1.4 


1.7 


10.5 
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PARANA! 




Figure 10. Nighttime wind roses at Paranal (1998-2006), and 
CAMC, TNG and NOT (1998-2004). 



Table 11. Nighttime yearly evolution of the wind direction fre- 
quencies at Paranal from 1998 to 2006 as [% per year]. 



N 


NE 


E 


SE 


S 


SW 


W 


NW 


-1.6 


+0.3 


+0.2 


+0.9 


+0.5 


+0.3 


< 0.1 


-0.3 



oscillation of the wind from north until 1994 (see Figure 
1 1 1 P . Data from 1998 show a clear trend as shown in Table 
[TT]that reports the computed trend of the yearly evolution of 
the frequencies in each direction. Wind from north is char- 
acterized by a decrease of the frequency of 1.6% per year 
marginally compensated by and increasing of wind coming 
from south-east (0.9% per year). We can conclude that pro- 
gressively wind coming from the sea is replaced by wind 
coming from the Atacama Fault. At the moment it is not 
clear if we are in presence of a wide-scale changing of the 
atmospheric conditions or if this effect is induced by local 
conditions. 

5.2 Wind speed 

The analysis of the wind speed at Paranal and ORM is car- 
ried out from the calculation of the time in which w sp is in 
fixed intervals established on the basis of the safety opera- 
tions at the observatories. We consider five main situations: 

• w S p < 3 m s _ : negligible wind speed, typically in this 
case the seeing increases (Paper II); 

• 3 ^ w sp < 12 m s _1 : the wind speed is in the safety 
range, telescopes observe without restriction in the pointing 
direction, seeing conditions are optimal; 

• 12 < w sp < 15 m s _1 : in this interval the telescopes can 
only point to objects in a direction 90° with respect to 
the actual Wdir)', 

• w S p > 15 m s _1 : at ORM the telescopes are closed for 
strong wind; 



Table 12. Nighttime wind speed statistics at Paranal (1998-2006) 
and CAMC, TNG and NOT (1998-2004). 



w sp range 
[m s" 1 ] 


Paranal 

[%] 


CAMC 

[%) 


TNG 

[%) 


NOT 
[%] 


w sp < 3 


22.1 


83.6 


30.2 


18.5 


3 ^ w sp < 12 


68.7 


16.4 


68.4 


70.2 


12 ^ w S p < 15 


5.9 


0.0 


1.1 


7.1 


w sp ^ 15 
■w sp ^ IS 


3.3 
< 0.5 


().() 
0.0 


0.3 
().() 


4.2 
1.2 



• w sp > 18 m s 1 : at Paranal the telescopes are closed 
for strong wind. 

Table [121 reports the percentages of time in which the wind 
speed is in the fixed intervals. Also in this case, because of 
the wind speed differences noticed between CAMC, TNG 
and NOT, none of them can be considered as fully represen- 
tative of the ORM, so they have to be taken into account 
together. For this reason Table [12] reports also the results 
for TNG and NOT obtained in Paper II. 
The site of CAMC has a predominance of w sp < 3 m s _1 
(83.6%), while Paranal has 22.1%, TNG the 30.2% and NOT 
the 18.5%. Paranal preserves good wind speed conditions 
in the ~70% of the time (same of TNG and NOT), while 
CAMC only the ~16%. CAMC never shows w sp > 12 m s" 1 , 
while Paranal has a percentage of 5.9 in the interval [12, 15[ 
m s _1 and 3.3% of w sp > 15 m s _1 while not-observing 
conditions due to strong wind occur in < 0.5% of the time. 
The TNG has ~1% of the time with w sp > 12 m s" 1 , while 
NOT has considerable higher values in the intervals [12, 15[ 
ms" 1 s" 1 (7.1%) and w ap > 15 m s _1 (4.2%). The highest 
wind speed measured at Paranal during the years 1998-2006 
is 27.4 m s" 1 in May 2000. 



6 CONCLUSIONS 

We have compared the astroclimatological conditions at the 
Paranal Observatory and El Roque de Los Muchachos Ob- 
servatory. Temperature, air pressure and relative humidity 
~20 years trends have been taken into account and signif- 
icant differences resulted form the analysis. In particular, 
CAMC denotes a warming of ~1.0°C/10yr, TNG data after 
1998 show a fiat trend, while the annual temperature trend 
at Paranal is characterized by a negative slope between 1993 
and 2000 (-0.12) while data in the period 2000-2006 have 
a positive slope (0.19). We have shown that using the pres- 
sure height scale both observatories are dominated by high 
pressure values and the computed pressure height scale at 
Paranal is 8334 m while at ORM it is 8325 m. 
Typically in cold seasons RH is lower than 60% at CAMC 
and 15% at Paranal. In warm seasons RH is lower than 
40% at CAMC and 20% at Paranal. A standard require- 
ment for the use of telescopes is a RH value < 80% or 
< 85%. The number of nights in which RH has been higher 
than the mentioned limits for more than 50% of the duration 
of the night denote a significant difference between Paranal 
and CAMC, the first appearing almost immune to high RH 
events. 

We have calculated the nighttime annual percentage of time 
in which the difference between the air temperature and the 
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Figure 11. Nighttime evolution of the wind direction through 
the years at Paranal (quadrant W-NE). 



Figure 12. Nighttime evolution of the wind direction through 
the years at Paranal (quadrant E-SW). 



dew point temperature is lower with respect to a variable 
upper limit X (X = 1 and 5°C). The analysis has shown 
better conditions at Paranal with respect to CAMC in win- 
ter, autumn and spring until 2000, while the two sites are 
almost equivalent in the respective summer. Since 2001 the 
two sites are becoming more similar. 

We have found a negative correlations between the Southern 
Oscillation Index and annual temperatures at Paranal (c.l. 
0.7), while no correlation has been found between SOI and 
yearly air pressure in the same site. 

Paranal is typically dominated by wind from north and 
north-east in nighttime. The wind speed is lower then 12 
m s^ 1 (safety range) in more than 80% of the time, while it 
is higher than 18 m s _1 less than the 0.5%. The mean wind 
speed is 7.0 m s _1 during the night. The same analysis at 
CAMC (see Paper II) shows that the mean wind speed is 
2.2 m s _1 , but significant differences we found in both mean 
direction ad wind speed for each different telescope at ORM. 
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